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Abstract. We express the generalized Cauchy-Stieltjes transforms (GCST) 
of some particular Beta distributions depending on a positive parameter A as 
A-powered Cauchy-Stieltjes transforms (CST) of some probability measures. 
The CST of the latter measures are shown to be the geometric mean of the CST 
of the Wigner law together with another one. Moreover, they are absolutely 
continuous and we derive their densities by proving that they are the so-called 
Markov transforms of compactly-supported probability distributions. Finally, 
a detailed analysis is performed on one of the symmetric Markov transforms 
which interpolates between the Wigner (A = oo) and the arcsine (A = 1) 
distributions. We first write down its moments through a terminating series 
3F2 and show that they are polynomials in the variable 1/A, however they 
are no longer positive integer-valued as for A = 1, = 00 (for instance A = 2) 
thereby no general combinatorial interpretation holds. Second, we compute 
the free cumulants in the case when A = 2 and explain how to proceed in 
the cases when A = 3, 4. Problems of finding a parallel to the representation 
theory of the infinite symmetric group and an interpolating convolution are 
discussed. 



Let A > and ii\ a probability measure (possibly depending on A) with finite 
all order moments. The generalized Cauchy-Stieltjes transform (GCST) of nx is 
defined by 



for non real complex z ([T7], [5], [II]). For A = 1, it reduces to the (ordinary) 
Cauchy-Stieltjes transform (CST) which has been of great importance during the 
last two decades for both probabilists and algebraists due the central role it plays 
in free probability and representation theories ([T], [S]). Moreover, CST were ex- 
tensively studied and they are well-understood in the sense that for instance, a 
complete characterization of those functions is known and one has a relatively easy 
inversion formula due to Stieltjes ([9]). However, their generalized versions are 
more hard to handle as one may realize from the complicated inversion formulas in 
|17j . [5], [13]. In this paper, we adress the problem of relating both generalized and 
ordinary transforms, that is, given /i^ check whether there exists a probability 
measure vx such that 
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1. Motivation 




(1) 
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for z in some suitable complex region where {z — x)^ is single- valued and is ana- 
lytic (but we prefer taking the principal determination) and characterize v\ in the 
affirmative case. Doing so will lead for instance to the invertibility of 

1 1 

for z belonging to some neighborhood of infinity ([2]) and as a by-product to a kind 
of A- free cumulants generating function for /i^, referring to the case A = 1 (il6j). 
However, this may not be always possible as we shall see, that is, v\ may not be 
a probability measure for some A. Moreover, when v\ is shown to be a probability 
measure, it is not easy to check whether it is absolutely continuous or not and more 
harder will be to write down its density when it is so. This was behind our willing 
to get an insight into the above problem since GCST naturally appear in the study 
of a particular class of probability measures ii\ {[B]). Those probability measures 
may be mapped via affine transformations into Beta distributions 

/3a,6(dx) = (2 - x)"(2 + a;)^l[_2,2]{a^)rfa; 

with parameters a,b > —1 depending on A, and their monic orthogonal polynomials, 
say {Pn)n, are the only ones that admit an ultraspherical type generating function: 



(2) V^P„^(a:)z"= ^, xesupp(MA), 



n>0 



ux{z){fx{z) - xY 



where ux, fx satisfy some technical conditions in a suitable open complex region 
near z = so that ([2]) makes sense. Due to the orthognality of P^, n> 0, one gets 
after integrating both sides in ^ with respect to fix that 

"^(^)" / 7TTT ^^^>^idx), 

Jr {f\[z)-xY 

for small z. Thus, if fx is invertible, the last equality is rewritten as 

u\[fx\z)] = (~^^^('^^) 

so that uxifx^) GCST of /ix- Fortunately, as one easily sees from the reminder 

below, f^^ may be identified with the CST (in some neighborhood of infinity) of a 
semi-circle law of mean and variance depending on A. 

2. Reminder and results 

Assume fix has zero mean and unit variance, then fix has a generating function 
for orthogonal polynomials of ultraspherical-type if and only if it is compactly- 
supported and belongs to one of the four families given by the following data ([6]): 

uxiz)^z^, fx{z) ^ ^^^^z + -,\> 0, 

2 z 
z A 1 

= l-(A/2)z2' /^(^) ^ 2 " + z ' ^ > 

z^ A^ 1 1 

ux{z) = , f^{z)^— -z± , = + -,A>l/2. 

l±Az/y2A^ ^' 2A-1 x/2A~T z ' 

We did not displayed the density of fix for sake of clarity and we will do it later. 
Nevertheless, it is worth recalling that the first family corresponds to the monic 
Gegenbauer polynomials (C;^)n, the second one corresponds to (C^~^)„ and is 
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related to the Poisson kernel, while the remaining families correspond to shifted 
Jacobi polynomials whose paremeters differ by 1. Now, the inverse in composition's 
sense of the CST Gm,a of a semi-circular law of mean m and variance cr, say Km. a, 
is given by ([4j) 

(3) Km.A^) := G-'Az) ^ ^ 



Then, the reader can check our previous claim about />,. Henceforth, we will only 
make use of G := Go,i since the elemantary identity holds 



Gra,a{,z) 



1 



-G 



The first main result may be stated for sufficiently large z as 



(4) 



1 



{z — x)^ 
where a(A) + 7(A) = 1 and 
iixiz) = z^ 
u\{z) = 



l/A 



1 



u\{z) = 



l±z' 



a{X) 


= 1, 




a{X) 


^ 1 - 


\,G{z) = 


a{X) 


= 1 - 





1 



1 



z±2 

for the four families respectively. It follows by the characterization of CST of 
probability measures that 

1 1 

(z ~ x)^ 

defines a CST of some probability measure I'x for A provided that < a(A) = 
1 — 7(A) < 1. Note that under this condition, the CST of iy\ is the geometric mean 
of G, G and that one discards the values A G]1/2, 1[ for the second family for which 
G is the CST of the arcsine distribution (^16j). 

The second main result states that under the same condition, i'x is an absolutely 
continuous probability measure and gives its density. This follows from the fact 
that I'x is the so-called Markov transform of some compactly-supported probability 
measure tx ([9J), that is 



(5) 



/ -^—I'xidx) = cxp- / log( 



z — x)Tx{dx), 



and from Cifarelli and Regazzini's results (|9| p. 51). For the first and the second 
families, vx is given by the Wigner distribution and a symmetric deformation of it 
respectively, while for the remaining ones, it is a non symmetric deformation of the 
Wigner distribution. Since the latter is a universal limiting object (representation 
theory of the infinite symmetric group, spectral theory of large random matrices, 
free probability theory), we give a particular interest in vx corresponding to the 
second family. We first express its moments by means of a terminating 3F2 series 
interpolating between the moments of the Wigner and the arcsine distributions 
(Catalan and shifted Catalan numbers, [5^ p. 64). Unfortunately, the moments are 
no longer positive integer-valued as it is shown for A = 2 therefore no general com- 
binatorial interpretation holds. Nevertheless, they are polynomials in the variable 
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1/A. Finally, we use ([4]) to compute the inverse of its CST and the free cumulants 
generating function in the case A = 2 which involve a weighted sum of the Catalan 
and the shifted Catalan numbers. For A = 3, 4, this is a more complicated task and 
by Galois theory, this is a limitation rather than a restriction since we are led to 
find a root of a polynomial equation of degree A. 

Remarks 2.1. 1/For the second family and for the discarded values of X such that 
the condition < a(A) = 1 — 7(A) < 1 



(z — a;)- 



-fj.x{dx) 



1/A 



does not define the CST of a probability distribution. 
2/ It follows from (QP and that 



1 



{z — x)- 



-li\{dx) — exp—X / \og{z — x)T\{dx) 



-v\{dx) 



Similar identities already showed up in relation to Bayesian statistics p. 59). 

Throughout the paper, computations are performed up to constants depending 
on A, which normalize the finite positive measures involved here to be probability 
measures. The paper is divided into five sections: the first four sections are devoted 
to the four families /ia while the last one is devoted to the particular interest we 
give in the probability measure vx corresponding to the second family. 



3. Markov transforms: symmetric measures 
3.1. First family. On the one hand, 

^2 \ 



yi,x{dx) cx 1 



2(1 + A) 



and its image of under the map x ^ ■\/(l + X)/2x has the density proportional to 
On the other hand, it is easy to see that 

. . 

so that 



^Ai-y-^-T--"o.V(TTI)7i' 



+ A 



It follows that 



{'i-x^)^~'^^dxoc= [G{z)y 



(z - xY 



-dx 
'2^ 



Using the fact that the Wigner distribution is the Markov transform of the arcsine 
distribution (see [9] p. 64), one finally gets 
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Proposition 3.1. For A > 

1-2 

(6) 



{z-x) 



'i — x^)^ ^^"^ dx cx 



J _2 Z — X In 

exp-A / log(z-a;) — dx. 
J-2 7rv4 — x-^ 



Remark 3.1 (Gauss hypergometric function). Recall that the Gauss hypergeomet- 
ric function of a complex variable z is defined by 

{a)n(h)nZ'' 



2Fi{a,b, c;z) = ^ 



n>0 



(c)„ n\ 



z < 1 



where c G M \ Z_ and (a)„ is the Pochhammer symbol defined by: 
(a)o — 1, (a)ri = (a + n — 1) • • • (a + l)a, n> 1. 
When c > b > 0, 2F1 admits the integral representation: 



iFi{a,b,c;z) 



T{c) 



r(c-6)r(6) Jo 

Then, it is an easy exercice to see that 



c-6-1 



dt. 



1 



-2 (2:-a;)' 



2NA-1/2 



1 



CX 



{z + 2y 



Fi ( A,A + i 2A+ 1; 



2 + 2 



for enough large z. To recover one uses the identity 
2F1 ( A,A + i 2A+ l;z ) cx ^ 



2A 



for z in the unit disc. 

3.2. Second family. The density of ^\ reads 



Hx{dx) « ( 1 - ^ 



2x A-3/2 



(a;)da;, A > 1/2, 



and we map it using x ^ Xjlx to 

(4-x2)^"'/'l[_2,2](x)d:E, A>l/2. 

Now, 



so tliat one gets 



-2 - a;)-- 
Now recall that, for z e C \ [-2, 2], 

7-2 ^ - a; 
and that G'^{z) + 1 = zG{z). Then 



da; z — a/z^ — 4 



27r 



iiA(G)(z). 



+ - 4' 



1 - G2(z) = 2 - zG(z) = 2- 



2z 

z + - 4 



= Vz^~4G(z). 
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Using (yL6j) 

TT J_2 z - a; V4 - a;2 - 4 2 

for suitable 2;, one gets 
Proposition 3.2. 

(7) f -^,{^~x')'-"' dx^^^=^ = G'-\z)G...Uz)■ 

Moreover, 

rxidx) =(l-\) ^^^l[-2,2]{dx) + ^^-l±^idx), 

and is not a probability measure unless A > 1. 

Next, for A > 1, 

niz) := Gi-VA(^)Gl/^„(z), z e C \ [-2,2], 

is the CST of a probability measure vx. This is readily checked using Lemma II. 
2.2 in |15j. More precisely, one has for > 



^Tg[G'-'^\z)Gii^,,Jz)] = (1 - ^) arg[G(z)] + j arg[Garcsi„(2)] e] - 7r,0[ 



so that T\ is of imaginary type (maps the upper half-plane into the lower half- 
plane), and 

lim lyTxiiy) = lim [^yG(^y)]l~l/^[^yGarcsi„(^y)]'/^ 

y — ^00 y — >oo 

For 1/2 < A < 1, one easily gets the inequality 

2arg[Garcsin(z)] < arg[TA(z)] < arg[Garcsin(2)] - arg[G(z)] 

and there is no guarantee to that vx is a probability distribution. In order to 
characterize ^a, A > 1, the Markov transform of tx, we will use results by Cifarelly 
and Regazzini p. 51). In fact, since tx is a compactly-supported probability 
measure, then vx is absolutely continuous with density proportional to 

(8) sin(7rFA(x)) exp — / log \x — u\Tx{du), 

where 

-!tFx{x) := taQ - co,x]) 

if a; < -2, 

7r/(2A) if a; = -2, 

(1 - l/A)[arcsin(a;/2) -|-7r/2] if a;e[-2,2[, ' 

TT if a; > 2 

Note that since ([8]) is valid when tx is the arcsine distribution and I'x is the Wigner 
distribution, one deduces that 

11 

exp— / logja: — uj d.ii. a;G[— 2.2], 

J -2 TT V4 - 1*2 
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does not depend on x (is constant). This striking result may be used to derive the 
density of vx in our case and in the forthcoming ones since t\ is a convex hnear 
combination of the arcsine distribution and a discrete probability measure. In the 
case in hand, easy computations yield 



Proposition 3.3. 

vxjdx) 
dx 



oc cos 



1 \ . X 
1 — — I arcsm — 



(4_^2)l/2A 1]-2.2[W,A>1. 



Note that A = 1 corresponds to the arcsine distribution while X = oo correponds 
to the Wigner distribution. Thus, v\ interpolates between them. The reader may 
wonder how to compute the normalizing constant or the moments of vx ■ This will 
be clear after dealing with the two remaining families ^x- 

Remark 3.2. One can derive ^ using the Gauss hypergeometric function. This 
time, 



1 



2N A-3/2 



1 



OC 



for enough large z, and one makes use of 
2F1 ( A - i A,2A- l;z I cx 



Fi ( A,A - i,2A- 1; 



1 



+ 2 



for z in the unit disc since 

4. Markov transforms: non symmetric measures 
4.1. Third family. The probability distribution ^x has the density 
V2A- Ix- / V2A-la;-l^ 



1 



2A 



1 



where 



Moreover 



X e 



1 - 2A 1 + 2A 



V2A- 1' V2A- 1 



2A 

, A > 1/2. 



A-3/2 



/a(^) = 



-z + 



2A- 1 V2A~ 



Thus, 



fx ^i^) — ^l/V2X^,X/V2X^(^) — ^ ^ 



A 

V2A- 1 
A 



G 



V2A- 1 
A 



V2A- Iz- 1 
A 



1 



V2A-1 



Now, the image of fix under the map x i-^- {\/2X — Ix — 1)/A transforms its density 
to 



1 



A-1/2 



On) 



A-3/2 



and one easily sees that 

'—U-x^y~"'{2-x)dxoc- 
xV^ ' ^ ' 1 



{z - xY 



G{z) 



ux{G){z). 
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Now note that (1-G(z))2 = l + G^{z)-2Giz) = (2- 2)G(z) which yields l-G(z) = 
\fz^^\fG{z) where the branch of the square root is taken so that 1— G is of positive 
imaginary type (i.e. maps the upper half plane into itself). As a result, 



Proposition 4.1. 



1 



-2 - ^) 
In this case t\ satisfies for suitable z 



G^{z) _ G^-^l^{z) 

-G{Z) ^ y/I^ ■ 



/ 



\og{z - x)Tx{dx) = - 1 - - log(G(z)) - - log 



2A 



2A 



z-2 



whence we deduce that 



which is a probability measure for all A > 1/2. Besides, the same arguments used 
before show that for suitable z 

Gi-V(^^)(^) 

(Z-2)1/2A ' 

is the CST of a probability distribution Px,X> 1/2 which is absolutely continuous 
of density given by 

Proposition 4.2. 

uxidx) 



dx 



(X sm 



(arcsm- + - 



1 



(2-a;)V2A 



l]_2,2[(a;),A> 1/2. 



4.2. Fourth family. The density of n\ is given by 



1 - 



\/2A - la; + 1 



2A 



A-3/2 



1 + 



where 



X e 



-1-2A -1 + 2A 



s/2\ - la; + 1 
2A 



A > 1/2. 



A-1/2 



V2A- 1' V2A- 1. 
The density of the image of fix under the map x 1— > (\/2A — la; + 1)/A reads 

(1 - 1)'"'^' (1 + 1)'"'^' l[_,,,](a;) cx (4 - x'f-'" (2 + a;) l^.,,,^{x). 



The same scheme used to deal with the third family gives 

Tx{dx) 



2Ay -K^/A^ 



'^[-2,2]{dx) + ^5_2{dx) 



and that the CST of vx is given by 



Gi-V(^^)(z) 

(2 + 2)V2A • 

However, the density of v\ is somewhat different from the one in the previous case: 



vxjdx) 
dx 



(X cos 



1 — ^1 arcsin — 
2A 



(2 + ,)i/2A ll-2,2[(a^).A>l/2. 
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Remarks 4.1. 1/For the four families, the measure t\ may be mapped to 
1 1 



-A[^i^i]{dx) + Mx5-i{dx) + Nx5i{dx) 



TT VI- x^ 

for some positive constants M\,N\. A more wider class of measures including the 
above one were considered in [12] . 

2/For Q<a<l,b— 1 — a, one can always define an operation 

where i>i, 1^2,1^ o-fe probability measures. However, every probability measure will be 
idompotent and the operation is not commutative unless a — b — 1/2. 

5. On the MOMENTS OF THE SECOND MARKOV TRANSFORM 

It is known that the Wigner distribution is a universal limiting distribution: it 
is the spectral distribution of large rescaled random matrices from the so called 
Wigner ensemble, the limiting distribution of the rescaled Plancherel transition of 
the growth process for Young diagrams ([10 ) and more generally of the rescaled 
transition measure of rectangular diagrams associated with roots of some adjacent 
orthogonal polynomials ( 11 )• It also plays a crucial role in free probability theory 
where it appears as the central limiting distribution of the sum of free random 
variables ([16], [1]). Note also that the standard arcsine distribution is the central 
limiting distribution for the t-convolution with t = 1/2 ([3], [5], [T3]) and is the 
limiting distribution of the so-called Shrinkage process f|10|V 

As a matter a fact, it is interesting to find a parallel to the above facts when the 
Wigner law is replaced by the Markov transform iy\ corresponding to the second 
family. More interesting is to define a convolution operation that interpolates the t- 
convolutions for t = 1/2 (X = 1) and t = 1 (A = 00, free convolution) having iy\ as a 
central limiting distribution ([7]). Since v\ is symmetric and compactly-supported, 
it is entirely determined by its even moments and we claim 

Proposition 5.1. The normalizing constant of v\ is given by 

-r(l-l/(2A))r(3/2-l/(2A)) 



r(2-i/A) 

and the even moments of v\ may be expressed as 



2«„ M.A-o2n P /-n, 1 - 1/(2A), 3/2 - 1/(2A) 



2 



x'''vx{dx)^2'^\F2[ '-'^ 2 i/A,l :1 



Moreover, the moments are polynomials in the variable (1/A). 
Proof, make the change de variable a; i— > 2 sin x in the integral 



x^" cos 



2 



1 — — arcsm 



A / 2 



to obtain 



22n+2-lA /"^'[sind^", 



1 

1 - - 



[A~ x^]-^/^^^Ux, n > 0, 

cosxl^-^/^dx, n > 0. 



Then, expand 



A 



[sin.x]2" = (1 - cos^x)" = ( J(-1)''[cos:e]"' 
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and use the formula (see [17], p. 177) 

(9) / coa\(p-q)x]\cosx]P+'^~^dx-~- ^ , 
with p + q = 2fc + 3— 1/ X,p — q = 1 — 1/A, to get 



TT r{p + q - 1) 



,p + q>l, 



s^" cos 



1 — — I arcsin — 



^ J ^4-^ Uy r(fc + 2 - 1/, 



-1/A) / 1 
- 1/A)fc! ^ 4 



This may be expressed through 3F2 hypergeometric series as follows: write the 
binomial coefficient as 

n\ ^ n\ ^ ..k i-n)k 
kj k\{n-ky. ^ ' kl 

and use the duplication formula to rewrite 

V^T (^2fc + 2 - 1^ ^ (^fc + 1 - ^) r (^fc + 3/2 - ^ 

It follows that 



x^" cos 



1 M • ^ 
1 — — I arcsm — 



[4 - xT'/^^dx . C.2-3F. 1 - 1/(2^)^3/2 - 1/(2A)^ 



To prove the last claim, let y 1/A and expand: 

(-f)..^(-f)('-'-f)-(-|) 



--) (y-2)(y-4)...(y-2fc + 2)(y-2fc), 



and similarly 



2 2 



(y-3)(y-5)---(y-2fc-l), 



r 

(2-y)fe = (-l)'^(2;-2)(2;-3).--(y-fc-l). 
The proof ends after forming the ratio 



(l-2;/2)fc(3/2-y/2)fc 



(2-y)A 



so that 



^"(^)-E(;)4 

k=0 



-\j iy~k~2)---{y-2k)iy~2k-l) 

,,_,{ y-k-2)---{y-2k){y-2k-l) 
k\ 



Remark 5.1. The exponential generating series 0/ (p„)„ is given by 

( \^ 4"-^ (j/-fc-2)...(y-2fc)(j/-2fc-l) 

2^P^^y>n\ ~ ^ (n-A:)!fc!' 

«>0 0<fc<ri ^ ' 



k\ 



^4, ^ - fc - 2) ■ • • (2/ - 2fc)(2/ - 2fc - 1) 

k>0 



k\ 



kl 



2'-ye'^,F,(^-y/^\(l-yy^^;Az) 
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5.1. Combinatorics. For A = 1, one recovers the moments of the arcsme distri- 
bution given by (use duplication formula) 



. (l/2)n 



When A = oo, one recovers the moments of the Wigner distribution: 



1 



n + 1\ n 



2n 



22„(l/2), 



(2)„ 



Both moments are integers and it is known that they count the shifted and the 
ordinary Catalan paths respectively (see |9j, p. 64). Unfortunately, TO2n 
positive integer-valued in general therefore no similar combinatorial interpretation 
can be given for all values of A. For instance, when A = 2, one gets 



k=Q 



(-1)* 



1 r(2fc + 3/2) _22„-i 



kj y/^ r(2fc + 2) 



fc=0 



4fc + 2\ 1 



2fc + 1/ 24fc 



by Gauss duplication formula so that the few first even moments are 3/2, 31/8, 187/16, 4859/128. 

Remark 5.2 (Non symmetric Markov transforms). The computations of the mo- 
ments of the non symmetric Markov transforms need more care. Consider for 
instance the one corresponding to the last family, then one has to compute after the 
change of variables a; > 2 sin x 



5Ji+l-l/(2A) 



smx cos 



-7r/2 



1 - 



1 

2A 



(1 + sinx)i/2A 



dx. 



To proceed, use the binomial Theorem to expand 

1 



(l-hsina;)i/2A 



/c>0 



and Fubini's Theorem to exchange both integral and sum signs. In fact. 



Tr/2 



smx cos 



Tr/2 



1 - 
7r/2 



1 

2A 



cosx > ' ' I sinxl'^dx < 
ki 

k>0 



1"^ cosx 
_./2 (l-|sinx|)i/2A 



dx 



cos X 



-dx 



[(l-sina;)i^i/2^]o/' 







(l-sina;)i/2A ~ 1 - 1/(2A) 

for A > 1/2. l^ow, it only remains to compute the integral 



< 00 



Tr/2 



Ism cos 



-7r/2 



1 



1 

2A 



cos xdx 



which is zero when n + k is odd and is evaluated similarly as the even moments of 
the above symmetric Markov transform when n -\- k is even using 
For the third family, make the substitution x )■ arcsin(a;/2) + 7r/2 in the integral 
(■2 



1 



1 
2A 



8 



to get 



2n+l-l/(2A)^_-^yi / [COSX]" 



' sm 



1 - 



2A 



1 



U/2A 



dx 



(1 -hcosa:)i/2A 



dx. 
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Then use the trigonometric formula 



sin a sin ^ = - [cos(a — 6) — cos(a + b)] 



to split the last integral to 
1 



/ [cosa;!"^"^- -r-7-T-dx — / [cos xl" cos 

Jo (l + cosa:)i/2A 2io ^ J 

Finally, use the binomial Theorem and 



2A 



(1 + cosxy/'^^ 



dx. 



5.2. Inverses of CST for A = 2, 3, 4. Let A = n > 1 be a positive integer. Recall 
that 

[G„(z)]" := [ / —L^v^{dx) = Un{G{z)), 
Uk [z-x) 

where as before G is the CST of the Wigner distribution. For the first family, 
Gn — G for all n so that Kn = K and K{z) = z + l/z for z in a neighborhood of 
zero. For the remaining families, the inverse of G„ is given by 

G-\z) X„(z) = if (fi-i(z")), K := G-\ 

for small |z|, subject to the condition 



Kn{z) = - + Rn{z) = - + 



TkZ 



k>0 

where R„ is an entire function known in free probability theory as the free cumulants 
generating function of z^„, {rk)k is the sequence of the free cumulants, rg, ri are the 
mean and the variance of Vn respectively. For the second family, one has to invert 

Z ^ Un{z) 



for small |z| in the lower half unit disc (image of G). We are thus led to find one 
root of the polynomial z" + z'^w — w for complex numbers z, w in a neighborhood 
of zero. This task is very complicated and even more, Galois theory tells us that 
the roots cannot be expressed by means of radicals unless n < 4. For n = 2, easy 
computations show that 



fe>0 ^ ' 



Thus one sees that tq = 0, ri = 3/2 which agrees with our above computations, 
and that 

■(l/2)fc , l(l/2)fe- 



r2k+i - (-1)* 



(1)a 



2 (2)s 



Thus, 



(-1)'=22'=+V2,+ 



22fe+i(l/2)^^22fe(l/2)fe 



Wk 



(2)a 



2k 



k + l\k 



For n = 4, one has a quadratic polynomial and setting v = z'^, one is led to 



— w = so that 
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The case z = 3 is more complicated and we supply one way to adress it: make the 
substitution ^ = Z — w/3 for suitable Z to get 

z"^ + WZ - w = Z"^ —Z + —w'^ - w. 

3 27 

The last polynomial has the same form as 

(a + bf - 3ab{a + 6) - (a^ + 6^) = 
which hints to look for a root of the form Z = a + b where 

2 



ab=—, a'^ + b'^ = w 



27 



and computations arc left to the curious reader. 

Remark 5.3. The above line of thinking remains valid for the non symmetric 
Markov transforms for which 

For A = 3, one already has the appropriate form of the polynomial and there is no 

need to make the above substitution. However, the case A = 4 needs more developed 
techniques since the polynomial is of degree four and is not quadratic. 
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